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ABSTRACT 


The design and performance analysis of a three-channel Robust 
Symmetrical Number System (RSNS) virtual spacing direction finding (DF) 
system was carried out in this thesis. The design was based on hardware used 
in the previous work in which a digital three-channel RSNS-based interferometer 
was built. In designing a RSNS virtual spacing DF system, it was found that the 
ratio between the modulus is an important parameter. The ratio affects the 
selection of the moduli set, the array size and the system’s field of view. Based 
on the hardware used, the system’s signal-to-noise ratio (SNR) was determined 
for a simulated emitter. Random Gaussian noise was injected into the system to 
determine the accuracy of mapping the actual angle-of-arrival (AOA) to the 
estimated AOA. Monte Carlo simulations were carried out to determine the 
probability of correctly estimating the AOA of the signal. Simulation results from 
the effect of spacing errors, phase errors, frequency error, moduli set change and 
SNR change was presented and analyzed. The RSNS-based DF system was 
found to be relatively insensitive to spacing error. A proposed two set of moduli 
were used to improve the accuracy and probability of correct estimation in both 
low SNR and high SNR conditions. 
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I. INTRODUCTION 


A. DIRECTION FINDING 

Radio direction finding is a class of direction finding (DF) by which the 
direction to a microwave signal is determined by means of a radio direction finder 
(RDF). A RDF is a receiving system that analyzes the information extracted from 
the passing electromagnetic radio wave to obtain direction of arrival (DOA) 
information [1], In military context, some of the more important uses of direction 
finding relating to electronic warfare applications include early warning threat 
detection, targeting and electronic intelligence [2], 

An RDF could be separated into three broad categories depending on the 
method by which it determined the angle of arrival (AOA) of the signal. The three 
categories are [2]: 

1. Amplitude Comparison 

2. Phase Delay 

3. Time Delay 

Amplitude Comparison uses either a direct or comparative method to obtain AOA 
information from the amplitude response of the antenna elements. 

Phase Delay and Time Delay systems extract the AOA information by a 
similar principle except that the former used phase information while the latter 
used the time of arrival information of the wave. Both systems required at least 
two separate antennas with pre-determined spacing. Figure 1 showed how a 
phase delay or time delay method would work. The plane wave would arrive at 
antenna element 1 earlier than antenna element 2. The antenna spacing, d, 
would cause a time delay or phase delay between the antenna elements. 
Assuming that the frequency of the wave and the antenna spacing are known, 
the AOA can be determined. 
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Incident Plane Wave 



Figure 1. Basic phase-delay / time-delay interferometry (After [1 ]). 


To increase resolution, a larger antenna spacing is usually used. 
However, a larger antenna spacing also introduces ambiguities. If the spacing is 
greater than a half-wavelength distance, more than one unique AOA exists. 

B. PREVIOUS WORK 

Through the years, research into using a Robust Symmetrical Number 
System (RSNS) for direction finding has been investigated in order to achieve 
high resolution with a small antenna baseline. It has been shown that a three- 
channel DF array based on RSNS theory and using a digital architecture can 
achieve a good resolution of 1.25 degrees RMS with a baseline of 66.2 cm [3]. 

C. PRINCIPAL CONTRIBUTIONS 

This thesis focused on the design and analysis of a three-channel virtual 
spacing DF system with two levels of processing. Design considerations for the 
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RSNS virtual spacing DF system was analyzed. It highlighted the importance of 
selecting the appropriate modulus ratio for sizing of the antenna array as well as 
field of view requirements. 

By analyzing the hardware used in the previous implementation of the 
RSNS-based DF system, the receiver’s noise figure was determined. From the 
noise figure, a link equation was used to calculate the theoretical signal-to-noise 
ratio (SNR) of a simulated emitter. Simulations were carried out using this link 
equation in determining the performance of a three-channel RSNS virtual 
spacing DF system. The effects of spacing errors, phase errors, frequency 
errors, moduli set change and SNR change were determined from the results of 
the simulation. Monte Carlo simulations were carried out to determine the 
probability of correctly estimating the AOA of the signal. Simulation results from 
the effect of spacing errors, phase errors, frequency error, moduli set change and 
SNR change was presented and analyzed. The RSNS-based DF system was 
found to be relatively insensitive to spacing error. A proposed two set of moduli 
were used to improve the accuracy and probability of correct estimation in both 
low SNR and high SNR conditions. 

D. THESIS OUTLINE 

Chapter II reviews the basic operating principles of RDF systems and 
provides the fundamental basis for applying the RSNS for direction finding. 

Chapter III introduces the RSNS theory and how it can be applied to 
phase interferometry. Previous experimental results are also presented. 

Chapter IV explains how the RSNS interferometry works with virtual 
resolution. 

Chapter V presents the considerations for designing a three-channel 
RSNS virtual spacing direction finding system. 

Chapter VI explains the models used in the simulation. 

Chapter VII presents the results and analysis of the simulations. 

Chapter VIII presents the conclusion. 
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II. PHASE DELAYED-BASED INTERFEROMETRY 


This chapter reviews the basic operating principle of phase delayed-based 
interferometry and how AOA ambiguities can arise. 

A. THEORY 

It is assumed that the DF array operates in the far field; hence the incident 
electromagnetic wave is planar. A one-channel linear interferometer is shown in 
Figure 2. The two antenna elements are space a distance d apart and the 
incident plane wave arrives at an angle 0 as shown. The angle of arrival can 
vary over a range of ±90°. 



Vout 

Figure 2. Two-element interferometry (From [4]). 
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The incident wave arrives at antenna element 1 first and after an 
additional distance of ds\n{0), it arrives at antenna element 2. A sinusoidal 
(time-harmonic) signal received by the antenna elements has the following form: 

Vj = V cos (co c t + y/j + (j>j) (2.1) 

where / denotes the element index number; V \s the signal amplitude response, 
co c is the signal frequency in radians/sec, y/j is the phase term due to the cables 


and (t>j is the phase term due to path delay relative to the reference antenna 

(antenna 1). The received signals are fed into an analog phase detector which 
consists of a mixer and a low-pass filter (LPF). After the phase detection mixing 
process, the output voltage has the following form: 

i/2 

V mix {</>J) = Vi{<faJ)V 2 {foJ) = —[cos{2a/ c t + J]y/ + Y J </>) + COS{Ay/ + A0)] (2.2) 

The term Ay/ results from the phase differential due to cable length differences. 
The term A <j> results from the phase differential due to the additional d sin(<9) 
distance traveled by the signal wave-front to reach antenna element 2 and it can 
be written as 


A(j> = ^ -(j >2 = kds'\r\{0) 


where k = 


2k 

T 


and X is the wavelength. 


(2.3) 


After the LPF removes the high frequency component of the signal, the 
output signal which remains is 


V* 

V ouM) = —COS(Ay/ + A(f) 


(2.4) 


We can measure Ay/ and eliminate it. After normalizing the voltage’s amplitude, 
the final output signal is 

Vout{</>) = cos{A</>) = cos{ kd s\r\{0)) (2.5) 


6 



B. AMBIGUITY 


As the phase differential term repeats itself with a period of 2 n, it can 
result in an ambiguity if distance between the antenna elements is greater than 

half a wavelength. To illustrate this point, we first consider a case where d = , 

which gives A<f> = ^sin(^). There will be a one-to-one mapping of the AOA to the 
phase difference as shown in Figure 3 and no ambiguity exists. 



Angle of Arrival (degrees) 


Figure 3. Phase differential vs AOA for d = j/^ (From [5]). 

Next, we consider a case where d = A, which gives A</> = 2^sin(^). For 
every phase difference measured, there will be two possible AOAs as shown in 
Figure 4. 
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Angle of Arrival (degrees) 


Figure 4. Phase differential vs AOA for d=X (From [6]). 

The accuracy of the AOA measurement can be stated as the change in 
phase difference for a given change in AOA. Flence from Equation (2.4), 

d{A<f>) = fa/QQs^) (2.6) 

86 

For any given frequency and AOA, Equation (2.6) shows that the accuracy of the 
measurement is directly proportional to the distance d. Flence, a longer baseline 
would give a more accurate measurement. 

To overcome the ambiguity coupled with the need for higher angular 
resolution, dual baselines are often used. The longer baseline is used to provide 
the high angular resolution required, while the shorter baseline is used to resolve 
any ambiguities [2], 

C. FOLDING WAVEFORMS 

The output voltage from any two antenna elements can be represented by 
Equation (2.5). For the cases of d = j/^ and d = A, the normalized voltages are 
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plotted in Figure 5. It shows that the both waveforms are symmetrical about the 
zero angle of arrival (broadside) and the visible region is from ±90 degrees. The 
number of folding periods, n, within the visible region is given by 



(2.7) 



-100 -80 -60 -40 -20 0 20 40 60 80 100 

Angle of Arrival (degrees) 



Figure 5. Output voltage vs AOA fore/ = ^ and d = A (From [3]). 

In a binary interferometer system where each successive baseline is 
changed by a factor of two, it is possible to resolve an ambiguity of the longer 
baseline and yet achieve higher resolution [2], A three-channel interferometer is 
shown in Figure 6. The successive baseline is changed by the factor of 2. 
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Y Y Y Y 

*- d 2 =A -* 


d,=24 


Figure 6. Three-channel interferometer with binary relationship 

between spacing (From [2]). 


The size of the physical antenna element array is often the limitation in 
determining the accuracy of the DF system. The next chapter introduces another 
method of increasing resolution without the requirement of a larger antenna 
element baseline. 
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III. ROBUST SYMMETRICAL NUMBER SYSTEM-BASED 

DIRECTION FINDING 


As it is advantageous to have a higher resolution without the need for a 
large antenna baseline, a signal preprocessing method utilizing the robust 
symmetrical number system (RSNS) was introduced to DF in [3-7], This chapter 
reviews the RSNS theory and how it is implemented in phase interferomentry. 

A. THE ROBUST SYMMETRICAL NUMBER SYSTEM 

The main idea behind RSNS-based direction finding is to represent a 
symmetrical waveform out of an antenna by a RSNS sequence. The RSNS is 
based on a sequence defined as 

x m =[0,1,2,...., m-1,m,m + 1,...., 2,1] (3.1) 

where x m is a row vector and m is the modulus and is a positive integer (m>0). 

In an A/-channel RSNS [8], where N>2, the basic sequence for the /th 
channel with modulus m is 


x mi = [%> 0l >■■■■, 0/v,1o > 1l>----’I A/’-■ 1 A/] ( 3 - 2 ) 

The sequence is repeated in both directions and it forms a sequence with a 
period of 

Prsns = 2Nm i ( 3 - 3 ) 

Considering a single channel, the discrete states of the RSNS can be expressed 
as 


n-Sj 

N 


9 = \ 


2Nmj + N- n + Sj -~\ 
N 


Sj < n< Nrrij + s,- +1 

Nrrij + Sj + 2 < n < 2Nrrij + s,- -1 


(3.4) 


where 
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g is the nth level of channel /, 
rrij is the /th channel modulus, 

Sj is the corresponding sequence shift S/ =0,1,2, , A/ —1 (mod N) and 

N > 2 is the number of sequences in the system. 

The values (s 1 ,s 2 ,....,s A/ ) must form a complete residue system modulo N 

and because of the relative property of the shifts, one of the shift values will be 
set equal to 0. 

An A/-channel RSNS is formed by selecting N moduli and N shift values 
Sj, (1</< A/). Since the fundamental period for channel i is 2Nirij , the RSNS 

vectors must be a mutilple of 2A/m,-. Hence, the fundamental period for the 
RSNS is 

pf RSNS = {2m 1 A/,2m 2 A/,...,2m /v A/} (3.5) 

where {a 1 ,a 2 ,....,a/ V } is the least common multiple of a^,a 2 ,....,a N . From Equation 
(3.5), 

pf RSNS =2N{m [ ,m 2 ,..., m N } (3.6) 


To illustrate the RSNS, consider a modulus set of [3 4] where 
m-| =3(5! =0) and m 2 = 4(s 2 =1). Table 1 shows the sequences illustrating the 
inherent integer Gray-code property. 


/ 

0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

m-i = 3 

0 

1 

1 

2 

2 

3 

3 

2 

2 

1 

1 

0 

0 

1 

1 

2 

2 

3 

3 

2 

II 

E 

0 

0 

1 

1 

2 

2 

3 

3 

4 

4 

3 

3 

2 

2 

1 

1 

0 

0 

1 

1 

n 





0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 



Table 1. RSNS sequence for m 1 =3(5! =0)and m 2 = 4(s 2 =1) 

(From [6]). 
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The system dynamic range, M, is the maximum number of distinct vectors 
without any ambiguity. From Table 1, the dynamic range for [3 4] is M= 15. It 
was shown that the shifts employed will affect the dynamic range sequence 
starting and ending points but it does not have any effect on the value of M [9]. 

The dynamic range of several moduli in combination has been formulated 
in closed form [10]. For two channels, there are three cases of two relatively 
prime moduli m 1 and m 2 that yield M. 


Case 1: For m 1 > 3 and m 2 = m 1 +1, 


M = 3(m-| + m 2 )~ 6 

(3.7) 

Case II: For m 1 > 5 and m 2 = m 1 +2, 


M = 3(m-| + m 2 )-7 

(3.8) 

Case III: For m 1 > 5 and m 2 = m 1 +C, and C > 3, 


M = 4 m 1 +2m 2 -2 

(3.9) 


For the example in Table 1, which satisfies Equation (3.7), 

M = 3(m 1 + m 2 )-6 = 3(3 + 4)-6 = 15 . 

For three channels, there are two cases of closed form solution that have 
been formulated [11, 12], 

• Case I: For m- [ =2 k -X m 2 = 2 k +1 for m 1 > 3, 

M = +-^/n 1 +7 (3.10) 

• Case II: For m^ = 2 m 2 +1, m 3 = 4 m 1 +1 , 

M = 6 mf + 21 /77-| +3 (3.11) 

As the closed form solution covers only limited cases, it has been shown that the 
most efficient method of determining the dynamic range is by using a search 
algorithm [12, 13]. 
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B. THE RSNS INTERFEROMETER DESIGN 


The design of an RSNS interferometer has been simplified into eight steps 
by Lee in his Master’s thesis [3]. The eight steps are: 

1. Determine N, the number of channels required. 

2. Identify an integer valued modulus for each channel {m^ ...m N ). 

3. Determine the system dynamic range, M, based on the 
chosen moduli. 

4. Determine the number of folds,/?/, for each modulus. 


5. 


6 . 


The number of folds is determined by 


n, = ——— = ——— (3.12) 

P RSNS 2Nm i 

Determine the required antenna spacing, dj, for each channel. 


, X MA 

d; = n. — = - 

1 1 2 ANmj 

Re-mapping the Field of View (FOV). 


(3.13) 


This is an optional step. The gain of an antenna reduces as the 
AOA increases from the broadside and to get a high signal-to- 
noise ratio (SNR) for digital processing, the FOV can be 
compressed into a suitable range. This increases resolution and 
decreases the antenna spacing required [4], By keeping kand \Z 0 ut 
constant in Equation (2.5), we get 

^sin(<9 max ) = c/'sin(^ nax ) (3.14) 

where <9 max is the maximum FOV and d is the antenna spacing of 
the original configuration. The quantities 0' max and d' are the 

maximum FOV and antenna spacing for the scaled configuration. 
The scale factor, £, is defined as 
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(3.15) 


^ s ^ n (^max) 

d si n (^max) 

Combining Equation (3.13) and Equation (3.15), the scaled antenna 
spacing can be determined. 

d max = d ii = 4 (3-16) 

7. Determine the thresholds. 

For a digital architecture, the thresholds to amplitude analyze the 
folding phase waveform into the RSNS. It has been shown [4] that 
for a normalized folding waveform, the threshold for 
the yth comparator of modulus channel m,- can be determine by 

( ■ o 

m, - j + — 

Vj m .= cos - —7t , 1 < y < m,-. (3.17) 

J ’ 1 nij 

v ) 

When the waveform exceeds a threshold, it will be encoded to the 
next integer value. Figure 7 shows how the folding waveform in 
Table 1 can be encoded by Equation (3.17). At the boxed in 
section of M, the symmetrical residues are s-, = 0 (m 1 = 3) and 
s 2 =3 (m 2 =4). 


15 




Figure 7. Mapping RSNS folding waveforms to threshold values (From [5]). 


8. Calculate the phase adjustment term, £, for each channel. 


From Figure 7, both moduli do not have symmetry about the center 
of the dynamic range at n= 7. A phase adjustment term,^ , is thus 
required to shift the center of the dynamic range to the signal AOA 
at the broadside. From Equation (2.5), the phase adjustment term 
is added: 


V out (6>) = cos(kds\n6> + C) (3.18) 

The relationship between the RSNS coding and the AOA has been 
obtained at this point. Using this relationship, a transfer function of 
the estimated AOA against the true AOA can be plotted. Previous 
studies have shown that the estimated AOA, 0' u , for the i/thbin is 


given as [4]: 


Q\, = sin 1 


2u + 1 

Mi; 


1 _ 

€ 


and the resolution, r u , for the i/th bin is given by: 


r u = sin 


-1 


2u-M + 2 

~M4 


-sin 


-1 




2 u-M 

Mi; 


(3.19) 


(3.20) 
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c. 


DESIGN EXAMPLE 


A three-channel RSNS DF system was designed and implemented by 
Lee in his master’s thesis [3]. Using the eight steps as previously described, his 
design is used as an example to explain how these steps are carried out in detail. 

1. Determine N, the number of channels required. 

It was decided a three channel DF system was to be designed. Hence, 
N=3 and a total of four antenna elements would be required. 

2. Identify an integer valued modulus for each channel (m^...m N ). 

The moduli set of [7 15 29] was chosen for high resolution as well as a 
relatively short antenna baseline. It is important to note that there will be 
other viable alternatives to the moduli set and the current set was chosen 
because it consists of pairwise relatively prime (PRP) numbers. 

3. Determine the system dynamic range, M, based on the 
chosen moduli. 

Using a MATLAB search program, the dynamic range for the moduli set 
was found to be M= 444. 


4. Determine the number of folds, /i f -, for each modulus. 

Applying Equation (3.12), the number of folds for each channel sequence 
was found to be, 




M 

2A/m 1 


444 

2(3)(7) 


= 10.57 


n 2 


M 

2Nm 2 


444 

2(3)(15) 


4.93 


"3 


M 

2A/m 1 


444 

2(3X29) 


2.55 
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5. Determine the required antenna spacing, dj, for each channel. 

At this step, it was decided that a signal frequency of 2.4 GHz be selected 
as there were readily available component hardware in the laboratory, 
which could be used for building the array. For a signal of 2.4 GHz, the 
wavelength, A, is 12.5 cm. Applying Equation (3.13), the distances 
between pairs of antenna elements are, 

d, = ^1 = 10.57 


A 712 5^ 

d 2 = n 2 - = 4.93 —- =30.83 cm 
2 2 J 

A 5^ 

c/o = /io — = 4.93 =15.95 cm. 

2 V 2 J 

6. Re-mapping the Field of View (FOV). 

The FOV was not scaled and hence the scaled factor £=1 (±90°). 

7. Determine the thresholds. 

Each channel would have a number of thresholds corresponding to the 
chosen modulus number. For a modulus of 7, there would be 7 
thresholds. Similarly, modulus of 15 and 29 would have 15 and 29 
thresholds respectively. The threshold levels for the chosen moduli [7 15 
29] are shown graphically in Figure 8 below using Equation (3.17). 


12 51 

— =66.07 cm 
2 
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- -05 - 








J 

11 

ji 

m 2 =15 


0.5 4 

o 4 

-05 4 

-1 

- 15-1 - 1 

m 3 = 29 

Figure 8. Threshold levels for moduli set [7 15 29]. 

8. Calculate the phase adjustment term, £, for each channel. 

The phase adjustment terms, Ci > were similarly found by a MATLAB 
search program. The phase adjustment terms are shown in Table 2. 



Phase Adjustment Term 

Channel ^X^ 

-1.4960 radians 

Channel 2,^ 2 

0.0698 radians 

Channel 3, C 3 

0.8666 radians 


Table 2. Phase adjustment term for moduli set [7 15 29] (From [3]). 
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D. 


SIMULATION RESULTS 


Using Equation (3.18), the folding waveforms for all the channels were plotted in 
Figure 9. The corresponding encoded waveforms were also plotted in Figure 10. 
By using Equation (3.19), the estimated AOA can be plotted against the true 
AOA as shown in Figure 11. 



Figure 9. Folding waveforms with phase adjustment. 



Angle of Arrival (degrees) 

Figure 10. Phase adjusted folding waveforms after encoding. 
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Transfer Function for Moduli [7 15 29] 



Figure 11. Transfer function for demonstration array using 
moduli [7 15 29] (From [3]). 

E. EXPERIMENTAL RESULTS 

A RSNS-based DF array was built and tested using same design as 
described in Section D. The experimental results are shown in Figure 12. The 
experimental results showed that although there are errors appearing 
intermittently throughout the AOA, the basic feature of the transfer function 
remains similar to the ideal transfer function as shown in Figure 11. The errors 
are likely due to antenna spacing error and phase errors from the demodulator 
boards [3]. 
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Figure 12. Experimental results for demonstration array 
using moduli [7 15 29] (From [3]). 
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IV. RSNS INTERFEROMETRY WITH VIRTUAL SPACING 


Virtual spacing for RSNS interferomery was first proposed by Chen for a 
two-channel system and it was later extended to an A/-channel system by Lee [3, 
6]. This chapter introduced how quadrature demodulation works followed by an 
explanation of how a two-channel interferometer works. Finally, a three-channel 
virtual spacing DF system is presented. 

A. QUADRATURE DEMODULATION 

The bandpass signal can be expressed by the in-phase component and 
the qaudrature component as 

x band (0 = x l cos {o} c t) - x Q sin {o) c t) (4.1) 

where x, is the in-phase component, 

x Q is the quadrature component, and 


co c = In-f c , f c is the carrier frequency. 



Figure 13. Quadrature type demodulation of two-channel interferometer 

(From [3]). 
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The output voltage from the antenna elements as shown in Figure 13 can 
be expressed as 

V(P n ,0 = Re(A(p n )e'W>) (4.2) 

where n = 1,2,3 is the antenna element index, p n is the position at antenna 
element n, t is the time. A(p n ) is defined as 


A(Pn) = \^P„)\e’* [P " ) =MP n )+jAo(Pn) (4-3) 

The output voltage can be mixed with a local oscillator signal to get the in-phase 
and quadrature terms 

V l {p n ,t) = V{p n ,t)cos(co c t) 

1 1 (4.4) 

= - A l{Pn) + -[Ai{p n )cos{2Q) c t)-AQ{p n )s\n{2G) c t)] 

V Q {p n ,t) = -V{p n ,t)s\n{(D c t) 

1 1 (4.5) 

= --[Ai{p n )s\n{2ojct) + AQ{p n )cos{2oj c t)] + -A Q {p n ) 

After low-pass filtering, the baseband terms that remain are 

V/iPn) = ^A;{p n ) = ^\A{p n )\cos{</>{p n )) (4.6) 

V Q {p n ) = ^AQ{p n ) = ^\A{p n )\s\n{<f>{p n )) (4.7) 

For a plane wave, let |A(p n )|=1. After normalizing Equations (4.6) and 

(4.7), the in-phase and quadrature terms of the received signals by antenna 
element n are 

Vin = Vi(Pn) = COS(</>{p n )) (4.8) 

V Qn = V Q(Pn) = s H</>{P n )) ( 4 -9) 

The phase differential term between antenna element two and three can 
be found by 

= </>(p 2 )-(/>{p 3 ) 

f i/ 

f Q3 


= arctan 


V V I2 ) 


-arctan 


V V I3 ) 


(4.10) 


= kd 2 sin 0 
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Similarly, the phase differential term can be found for any other antenna 
pair. The in-phase and quadrature outputs for any channel using antenna 
element ras the reference, can be defined as 

V/ N (0) = cos{kd N s\n0) = cos{A<f> rn ) 

= cos(^(p n ))cos(^(p r )) + sin(^(p n ))sin(^(p r )) (4.11) 

= V ln V lr + V Qn V Qr 

V QN {0) = s\n{kd N s\nO) = sin(A^ m ) 

= sin(^(p n ))cos(^(p r ))-cos(^(p n ))sin(^(p r )) (4.12) 

= V ln V lr - V Qn V Qr 

The phase adjustment term, is required for channel /'so as to map the center 

of the dynamic range to the broadside. From Equation (3.18), the voltage from 
channel / can be expressed as 

Vj(0) = cos(/c(d / sin(6 , ) + 4' / ) (4.13) 

Rearranging Equation (4.13), 

V N (0) = cos(/ccy / sin(6 , ))cos(^ / )-sin(/rc/ / sin(6 , ))sin(^' / ) (4.14) 

Vi (0) = V n {0)cos(Ci)- V Qi (0)sin(d) (4.15) 

B. TWO-CHANNEL VIRTUAL SPACING INTERFEROMETRY 

For a two-channel RSNS interferometer as shown in Figure 13, the moduli 
pair m 1 and m 2 is chosen. The dynamic range of the pair is M. By combining the 
terms as given in Equation (4.12) and Equation (4.13), the in-phase and 


qradrature term of the virtual spacing can be obtained: 

V'n = V/ 1 V /2 - Vq^Vq 2 = cos(/c(cy 1 + d 2 )s\n0) = cos(kd\ s\n0) (4.16) 

Vq-\ = Vq 1 Vj2 + V/-\Vq 2 = sin(/c(c/ 1 + d 2 )s\n0) = s\n{kd\ sin<9) (4.17) 

V/ 2 = V/ 1 Vj 2 + Vq^V Q2 =cos(/c(c / 1 -d 2 )s\r\0) = cos{kd' 2 s\r\0) (4.18) 

Vq 2 = Vq^Vi 2 - V n V Q2 = sin(/c(c/ 1 -d 2 )s\n0) = s\n{kd 2 sin<9) (4.19) 

Now suppose that the chosen moduli set [m 1 m 2 ] is given by 

m 1 = hm 2 (4.20) 
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(4.24) 

(4.25) 


The resolution at broadside is given by Equation (3.20). For AOA = 0° at 
broadside, the bin number, u, corresponds to 90. 


26 



r u= 90= sin 1 


= sin "* 


2u-M + 2 


MZ j 

2(90)-35 + 2^ 
35(1) 


-sin 


2u-M 

v MZ 


-sin "I 


2(90)-35 
v 35(1) j 


= 3.3° 


(4.26) 


From Equation (4.20), 





r i4) 

m' 2 = m\ 

ll -h) 

= m\ 

i-i 


3 m\ 


(4.27) 


There could be several modular pairs which satisfy Equation (4.21) and Table 3 
shows some of the possible pairs and their corresponding dynamic ranges. 


m\ 

7 

8 

9 

10 

m '2 

21 

24 

27 

30 

M' 

45 

51 

57 

63 


Table 3. Moduli sets and dynamic range for h = (From [6]). 


Suppose that the chosen virtual moduli set is [m\ m' 2 ] = [10 30] with 


M' = 63. The virtual scale factor can be found from Equation (4.18) or Equation 
(4.19), 


dMtl4 = (6.8 + 3.4) = ^§f 


4(2)(10) 


£'=1.04 

The resolution of the virtual moduli set is 


r u= 90 = s| n 


1 


2u-M' + 2 


= sin 


= 1.8 C 


-1 


-sin 


-1 


M'Z' 

2(90)-63+ 2^j 
63(1.04) y 


2u-M 


C/ii \ 


M'£ 


-sin ! 


2(90)-63 
63(1.04) 


(4.28) 


(4.29) 
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The virtual moduli set improved upon the resolution of the original array by 
almost 100%. In theory, super resolution could be achieved with a large number 
virtual moduli set [m\ m' 2 \. However, the larger the m\ or m' 2 , there would be a 

corresponding increase in the number of thresholds. As shown in Figure 8, these 
levels of threshold are closer together at the peaks and troughs of the folding 
waveform. Having the threshold levels closer can lead to quantization error in 
the presence of noise if the signal-to-noise ratio is low. It has been shown by 
Chen that in a low SNR environment, every moduli set would have regions where 
the transfer function “fades,” i.e., where gross errors tend to occur [6]. 

C. THREE-CHANNEL VIRTUAL SPACING INTERFEROMETRY 


The virtual RSNS DF can be similarly extended to a A/-channel system [3]. 
Let us consider a case of three-channel virtual DF array where the moduli set are 
related by 

m 1 = hm 2 = h 2 m 3 (4.30) 

The ratio h must satisfy the condition for three integer values of m,- to 


exist. Similar to the two-channel system, the virtual spacing may be extended to 
a three-channel system by defining d 3 as 


d' 3 = d^-d 2 -d 3 = 


AM 


4 N 


J_1_1_ 

v m i m 2 m 3 j 


4 = 


AM 


4 N 


1-ft-ft 

m l 


.2 ^ 


AM’ 

^ ~ 4Nm' 3 * (4.31) 


Rearranging Equation (4.18) with Equation (4.19) and Equation (4.30), we get 


1 + /7 

1-ft 


m 2 

m \ 


(4.32) 


1 + /7 


nj3_ 

m\ 


(4.33) 


V1 — ft — // y 

Depending on the ratio h, it is possible to select integer value of [m\ m 2 m 3 ] so 


as to form other virtual moduli sets. A simple algorithm can be implemented to 
post-process the AOA information using the virtual moduli set to increase the 
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resolution of the DF array. A simple algorithm flow with two levels of virtual 
resolution processing is shown in Figure 14. 


• Estimate AOA (eAOAl) from moduli set 

• Estimate virtual AOA (eAOA2 and eAOA3) for 
2 set of virtual moduli where resolution of 
eAOA2>resolution of eAOA3 



Figure 14. Algorithm for two levels of virtual resolution processing (After [3]). 
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V. DESIGN CONSIDERATIONS FOR THREE-CHANNEL RSNS 
VIRTUAL SPACING DIRECTION FINDING SYSTEM 


This chapter presents some of the considerations in designing a three- 
channel RSNS virtual spacing DF system. 

A. LIMITED REAL MODULI SET 


Let [m 1 m 2 m 3 ] be defined as the real moduli set where m 1 < m 2 < m 3 . In 


order for virtual resolution to work, the real moduli set must satisfy Equation 
(4.27). Let the ratio, R, be defined as 


As mentioned in the previous chapter, h<J/^ is required for three distinct 


positive integer values of the moduli set to exist. From Equation (5.1), R> 2 
follows. 


Based on these conditions, some of the possible real moduli sets which 
can be selected for a three-channel RSNS DF system are shown in Table 4 
below. 
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[m 1 m 2 m 3 ] 

M 

CM 

II 

DC 

\OJ 

II 

-c 

[2 4 81 

29 

[3 6 121 

41 

[4 8 161 

53 

[5 10 201 

65 

h = y 2 , R = 3 

[2 6 181 

59 

[3 9 27] 

86 

[4 12 361 

113 

[5 15 451 

140 

'Cj- 

II 

DC 

II 

-c 

[2 8 321 

101 

[3 12 48] 

149 

[4 16 641 

197 

[5 20 80] 

245 


Table 4. Real moduli set [m 1 m 2 m 3 ] and dynamic range for 

several values of h. 


B. LIMITED VIRTUAL MODULI SET 

Let [m\ m 2 m 3 ] be defined as the virtual moduli set where m\< m' 2 < m 3 . 
For the virtual moduli set to have positive integer value, it needs to satisfy 
Equation (4.32) and Equation (4.33). Let the ratios, R\ and R 2 , be defined as 


R^ 


m 2 

m\ 


R. 


2 - 


m\ 


1 + h' 

1 -hj 
1 + h 




(5.2) 

(5.3) 


For the virtual moduli set to have positive integer value, m\ must satisfy both 
Equation (5.2) and Equation (5.3). Hence m\ has to be the multiple of lowest 
common denominator of R\ and R 2 . For /?=%, R\ =3 and R 2 = 6. Since the 
lowest common denominator of R^ and R 2 is 1, there is no restriction on the 


value of m\. For h=V,, R^= 2 and As the lowest common 
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denominator of ft, and ft 2 is 5, m\ can be any multiple of 5. Table 5 lists some 
of possible the virtual moduli sets for selected ratios. 



[m\ m’ 2 m 3 ] 

M' 

h ~Yi 

[2 6 12] 

41 

[3 9 18] 

59 

ft, = 2, R 2 = 6 

[4 12 241 

77 

[9 27 54] 

167 

h= y$ 

[51012] 

89 

[10 20 241 

179 

R|=2, n 2 = 1 % 

[15 30 361 

269 

[20 40 48] 

359 

h = /4 

[33 55 60] 

853 

[66 110 1201 

1708 


[99 165 1801 

2563 

[132 220 2401 

3418 


Table 5. Virtual moduli sets [m\ m 2 m' 3 ] and dynamic range for 

several values of h. 


C. SENSITIVITY OF SIZE OF ANTENNA 

The size of an antenna array is determined by Equation (3.16). If the 
design frequency f is fixed at 2.4 GHz, the antenna size is only dependent on 
the selection of modulus, m,-, dynamic range M and scale factor Since the 

dynamic range M is determined by modulus my, the antenna size is only 
dependent on two parameters, my and 

The antenna spacing of the moduli set is listed in Table 6. It shows that 
the antenna array spacing is relatively insensitive to the selection of the real 
moduli set and more dependent on the ratio h. 
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[m 1 nn 2 m 3 \ 

d\ (cm) 

d 2 (cm) 

d 3 (cm) 


[2 4 8] 

15.1 

7.6 

3.8 


14.2 

7.1 

3.6 

£ = 1 


13.8 

6.9 

3.5 


13.5 

6.8 

3.4 

h= y$ 

mmEU 

30.7 

10.2 

3.4 


29.9 

10.0 

3.3 

£ = i 

[4 12 36] 

29.4 

9.8 

3.3 


29.2 

9.7 

3.2 

h =/ 4 


52.6 

13.2 

3.3 ! 

[3 12 48] 

51.7 

12.9 

3.2 

f=i 


51.3 

12.8 

3.2 j 

[5 20 80] 

51.0 

12.8 

3.2 


Table 6. Antenna spacing for [m 1 m 2 m 3 ] 
(scale factor £ = 1 and frequency at 2.4 GHz). 


Changing the scale factor to £ = 1.154 reduces the FOV to ±60° but it 
increases the array spacing slightly as shown in Table 7. 



[m 1 nn 2 m 3 \ 

di (cm) 

d 2 (cm) 

d 3 (cm) 

h 

[2 4 81 

17.4 

8.7 

4.4 


16.4 

8.2 

4.1 

£ = 1.154 

■Eiaiari 

15.9 

8.0 

4.0 

[5 10 201 

15.6 

7.8 

3.9 I 

h= y 

[2 6 181 

35.5 

11.8 

3.9 

ilEEESH 

34.5 

11.5 

3.8 

£ = 1.154 

[4 12 361 

34.0 

11.3 

3.8 

[5 15 451 

33.7 

11.2 

3.7 

h ='A 


60.7 

15.2 

3.8 

[3 12 481 

59.7 

14.9 

3.7 

£ = 1.154 

[4 16 641 

59.2 

14.8 

3.7 j 

■bW»g;!»W 

58.9 

14.7 

3.7 | 


Table 7. Antenna spacing for [m 1 m 2 m 3 ] 
(scale factor £ = 1.154 and frequency at 2.4 GHz).. 
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D. FIELD OF VIEW LIMITATION 

The FOV of the real moduli set is determined by the Equation (3.15) while 
the FOV of the virtual moduli is determined by Equation (4.21), Equation (4.22) or 
Equation (4.31). The critical parameter in determining the FOV is the ratio h. 
Table 8 and Table 9 list the FOVs of selected combinations of real and virtual 
moduli sets for £ = 1 and £ = 1.154. 



[m 1 m 2 m 3 ] 

[m\ m' 2 m' 3 ] 

FOV{ deg) 

FOV'(deg) 

\\ 

II II 

-C 



±90 

±69.9 


[9 27 541 

±90 

±64.8 

[4 8 161 

[4 12 241 

±90 

±75.6 



±90 

±69.0 

[5 10 20] 

[4 12 241 

±90 

±80.8 


[9 27 54] 

±90 

±72.1 

\co 

II II 

-C 

[2 6 18] 

[51012] 

±90 

±26.9 


mEESEm 

±90 

±27.1 

■EEESH 

H3E1EH 

±90 

±27.8 


[10 20 24] 

±90 

±27.9 

IEKEE3I 


±90 

±28.2 


WEHI 

±90 

±28.4 

\^t 

II II 

-C 

[2 8 32] 

[33 55 601 

±90 

±24.2 


[66 110 120] 

±90 

±24.2 

■EISEn 

■EEEHEm 

±90 

±24.6 


[66 110 1201 

±90 

±24.6 

[4 16 64] 

[33 55 60] 

±90 

±24.8 I 


[66 110 1201 

±90 

±24.9 


Table 8. Field of view for selected moduli set, £ = 1. 
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[m 1 m 2 m 3 ] 

[m\ m' 2 m 3 ] 

FOV (deg) 

FOV'(deg) 

h ~/4 

£ = 1.154 

[3 6 121 

[4 12 24] 

±60 

±54.4 



±60 

±51.6 1 

[4 8 161 

[4 12 241 

±60 

±57.0 


[9 27 54] 

±60 

±54.0 



±60 

±58.8 


[9 27 541 

±60 

±55.5 

h = /i 

£ = 1.154 

[2 6 181 

[5 10 121 

±60 

±23.1 


[10 20 241 

±60 

±23.2 



±60 

±23.8 


[10 20 24] 

±60 

±23.9 

[4 12 361 

[5 10 121 

±60 

±24.2 


MimKM 

±60 

±24.3 

4 = 1.154 


mmmm 

±60 

±20.8 


[66 110 120] 

±60 

±20.8 

[3 12 481 

[33 55 60] 

±60 

±21.1 i 



±60 

±21.2 

[4 16 641 

[33 55 601 

±60 

±21.3 


■laaroiKiiB 

±60 

±21.3 


Table 9. Field of view for selected moduli set, £ = 1.154. 


The DF system’s FOV is determined primarily by h and is relatively 
insensitive to the selection of the moduli set. As the ratio h decreases, the 
system’s FOV reduces accordingly. For a given operational requirement of an 
instantaneous FOV, the choice of the ratio h could be limited. 

E. SELECTION OF MODULUS RATIO 

Let the phase differences between every combination of element of a 
three-channel interferometer be defined as indicated in Figure 15. An ambiguity 
diagram can be used to determine which h ratio should be chosen. An ambiguity 
diagram is a phase-plane plot comparing the phase difference of an 
interferometer [14]. The closer the lines on the ambiguity diagrams, the poorer 
the interferometer will perform in determining the actual phase difference. To 
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further explain how the ambiguity diagram works, consider a simple three 
element interferometer with d^=3^j/Qand d 2 = The phase-plane plot of 

the corresponding $ and ^ 2 is given in Figure 16. 


4 3 2 1 

V V V V 


<k — > 

--► 

- <t\ 


Figure 15. Phase difference for a three-channel interferometer. 


6 n 


4 K 

fa 

2 n 


0 

0 2 n Ak 

fal 



Figure 16. Phase-plane plot for = 3 () ancl ^ 2=2 (/^) ( From t 14 ]). 
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Since $ and can be known only to within 2 n radian, Figure 16 can be 
normalized to 2 n as shown in Figure 17. The ambiguity diagrams for selected 

and h = are given in Figures 18 to 20. These figures 

showed that the higher h ratio gives the best performance in resolving ambiguity 
between the phases as the lines are further apart. 


ratios of h = y , h = 



Figure 17. Ambiguity diagram for = 3 (^ 2 ) ancl ^ 2=2 (%) ( From t 14 ]). 

In designing a RSNS virtual spacing DF system, it is important to select 
the highest possible h ratio. For any given h ratio and frequency, the size of the 
antenna array as well as the FOV of the system would not change significantly 
once the ratio has been selected The highest possible h ratio gives the best 
performance in resolving phase ambiguity. 
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Figure 18. Ambiguity diagram for h = ^. 
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Figure 19. Ambiguity diagram for h = 
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Figure 20. Ambiguity diagram for h = . 
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VI. MODELING AND SIMULATION 


To understand and predict the performance of a three-channel RSNS 
virtual spacing DF system, simulations in MATLAB were carried out. This 
chapter introduces all the models that were used in the MATLAB simulations. 

A. NOISE MODEL 

It is assumed that the noise in the system is random in nature and 
Gaussian distributed. The bandpass Gaussian random noise signal can be 
represented by 

W( t) = Re (| A(t)\ e j(p{t) e j( ° ct ) (6.1) 

where 

\A(t)\ is random process with Rayleigh first-order probability density 
function; and 

cp(t)- random process with uniform first-order probability over range 

( 0,2 n). 

Using Euler’s formula, Equation (6.1) can be written as: 

W(t) = W,cos(a> c t) - W Q (t)s\r\(a> c t) (6.2) 

where 

Wi(t) = \A(t)\cos(<p(t)) (6.3) 

W Q (t) = \A(t)\s\n(cp(t)) (6.4) 

From these equations and combining Equation (4.4) and Equation (4.5), 
the Gaussian noise model can be represented by 

V, N0 iJp m t) = V(p m t) + W l (t) (6.5) 

V Q Nni JPrvt) = V(p n J)+W Q (t) ( 6 . 6 ) 
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B. MONTE CARLO SIMULATION 


Monte Carlo methods are commonly used in simulating physical or 
mathematical systems when a closed form statistical solution is too complex to 
resolve. It uses the Monte Carlo algorithm which can be easily implemented with 
a computer to model random processes. The Monte Carlo simulation is carried 
out to determine the performance of the three-channel RSNS virtual spacing DF 
system with random white Gaussian noise. First, a review of the Monte Carlo 
algorithm is given below [15]. 

Let a random variable X consist of both the signal voltage V together with 
random noise W, centered on zero with variance of one. 

X = V + W (6.7) 

The probability of X>z for any given V where z is constant is 


P(z) = P^X>^ j (6.8) 

By generating a series of statistically independent and identically 
distributed Gaussian random variable w(n), where n=1 ,2,...,N times, Equation 
(6.8) can be estimated. A new set of random variables, x(n), can be found by 
adding a constant vto the w(n). 

x{n) = v+w{n) (6.9) 

Let a new random variable y(n) be defined as 

II, if x(n) > z 


( 6 . 10 ) 


The estimation of Equation (6.8) is 


P{z) = — Y y{ri) 

N „.i 

The mean, £(•), and variance, a 2 , of P(z) is 


( 6 . 11 ) 


, „ . i 'v -I 

E ( P < Z >) = 771 E (yW) = 77 WP < Z > = P < z > (6.12) 

v n =1 
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<t? (z) =£(p(z)-P(z))=£(p 2 (z))-P 2 (z) = 1p(z)(1-P(z)) 

Let the required outcome probability of the estimation, a , be defined by 

a < p||P(z) - P(z)| < 

where the relative error s is defined as, 

P(z)-P(z) 

£ =- 

P(Z) 

Using Chebyshev’s inequality, 

a? 

P( P(z) - P(z) < £) > 1- > a 

' ' £* 

Using Equation (6.13), 

lp(z)(1-P(z)) 

1--“-2-- a 

N> p (z)(1- p (z)) 

(1 -a)£ 2 


(6.13) 

(6.14) 

(6.15) 

(6.16) 

(6.17) 


Equation (6.17) determines the number of times the simulation must be run to 
satisfy Equation (6.14) for any given a and.?. The Monte Carlo algorithm is 
summarized by the flowchart in Figure 21. 
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Figure 21. Monte Carlo algorithm. 

C. HARDWARE SETUP 

The simulations were based on the hardware used by Lee in his three- 
channel DF antenna array [3], Figure 22 shows the block diagram of all the 
hardware used in the array. 
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Figure 22. Hardware block diagram of a three-channel DF antenna array. 


To determine the SNR of the I and Q data after the demodulator board, a 
calculation of the noise figure was done for one receiver channel. The 
specifications of the hardware used are given in Figure 23. 




Figure 23. Specifications for antenna array’s noise figure. 
The receiver’s noise figure as shown in the dashed box is given by 


c F 2 -1 
= F\ + —— + 

1 G| 


*3-1, F 4 -1 
G\G 2 G1G2G3 


(6.18) 
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where 


F\ = 


^ = cable loss from antenna to LNA = 0.5 dB, 


F 2 = noise figure of LNA = 1.7 dB, 


G 2 =gain of LNA = 23 dB, 

F 3 = y G ^ = cable loss from LNA to AD 8347 demodulator = 0.65 dB, 
F 4 = noise figure of AD 8347 demodulator =11.7 dB, 


G 4 = total gain of the demodulator baseband and RF amplifier. 
From Equation (6.18), 


14.79-1 


c 10 1.48-1 1.16-1 

Fx — 1.12 h -1-1- 

e 0.89 0.89x199.5 0.89x199.5x0.86 

F e = 1-77 

T e =(F e -l)r a =(1.77-1)290 = 223 K 
T s =T e + T a = 223 + 300 = 523 K 


(6.19) 

(6.20) 
(6.21) 
(6.22) 


The noise power is 


N 0 = k B T s B 


(6.23) 


= 1.38 x 10 -23 x 523 x 50 x 10 6 

= 3.61x10“ 13 
= -124.4 dB 

where B is the instantaneous receiver bandwidth (50 MHz). The minimum RF 
power that can be input to the AD 8347 depends on the sensitivity of the 
demodulator. From the specification, it ranges from -65 dBm to 0 dBm. 
Therefore the minimum power that the antenna dipole can receive is 

Free,mm =-65 + 0.65-23 + 0.5-6.55 =-93.4 dBm (6.24) 

The minimum SNR output from the AD8347 demodulator board is 

SNR mln = -93.4 - (-124.4) = 31.0 dB (6.25) 

Assuming that atmospheric absorption is negligible, for any isotropic 
emitter transmitting an average power of P t at distance R, the SNR can be 
determined by 
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SNR = 


4nR 2 N n 


(6.26) 


Based on the calculated SNR, the variance of the noise is determined from 


SNR = 


y?_ 

2a 2 


(6.27) 


where V) is the channel output voltage defined in Equation (4.16). For a one 
watt emitter, the maximum detection range of the current array is 


SNR m\n = 


4^Rmax 


1259 = 


4;rx R max 2 x3.61x10 13 


^max _ 13.2 km 


(6.28) 

(6.29) 

(6.30) 


D. ANTENNA PATTERN 

The antenna array used by Lee utilizes vertical printed circuit dipoles and 
was designed for a nominal frequency of 2.4 GHz. It has been tested in the 
anechoic chamber by Burgstaller to determine the dipole antenna pattern [16]. 
Figure 24 shows the measured element patterns for a dipole antenna at 2.4 GHz. 
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Figure 24. Dipole element pattern measured at 2.4 GHz (From [13]). 


The gain of the dipole element was measured to be 6.55 dB or 4.52 at 
broadside. The array estimates the AOA in the horizontal plane (H-plane). For 
the simulation, a simple cosine model was used to replicate the antenna pattern 
for various AOA 

G a = 4.52xcos(<9) (6.31) 

All the tools used in the simulation have been presented and the next 
chapter presents the results from simulations that were carried out. 
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VII. RESULTS AND ANALYSIS 


With all the tools of the simulation in place, a sensitivity analysis was 
carried out to find out how some of these errors might affect the performance of 
the RSNS-based DF system. 

A. BASELINE DESIGN PARAMETER 

A three-channel RSNS virtual spacing DF design was chosen as a 
baseline design before any sensitivity analysis could be carried out. Using the 
steps described in Chapter III, the design parameters for a three-channel RSNS 
virtual DF system were chosen. It was decided that the design would be based 
on the hardware shown in Figure 22; hence the design frequency of 2.4 GHz was 
selected. 

1. Determine N, the number of channels required. 

A three-channel DF system was to be designed. Hence, N=3 and a 
total of four antenna elements would be required. 

2. Identify an integer valued modulus for each channel (mi, m 2 
m 3 ). 

Before identifying the moduli sets, the ratio h must first be selected. 
The three-channel RSNS virtual spacing DF system in Lee’s thesis 
was based on h = Y 2 - The size of the antenna array built was 

about 14 cm. Experiment results shown that a number of gross 
errors appeared throughout the transfer function. Although the 
sources of errors were not fully investigated by Lee, it was 
suspected that the errors were due to mutual couplings as the 
antenna spacing between the antenna elements was small. It was 
decided that h = would be appropriate to reduce the effects from 
mutual couplings as well as to obtain a reasonable system’s FOV. 
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Based on h = j^, the real moduli set of [2 6 18] was selected so as 

to reduce quantization error when noise is present. The virtual 
moduli sets of [5 10 12] and [15 30 36] were selected to enable the 
DF system to have sufficient resolution. 

3. Determine the system dynamic range, M, based on the 
chosen moduli. 

From Table 4 and Table 5, the dynamic ranges for the all the 
moduli sets are listed in Table 10. 


Moduli set 

Dynamic Range 

[2 6 18] 

MReal = 89 

[51012] 

M'virtuaH = 89 

[15 30 36] 

M'virtual2 = 269 


Table 10. Dynamic ranges for baseline design. 


4. Determine the numbers of folds,/?/, for each modulus. 

The numbers of folds are given in Table 11. 
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Moduli Set 

Number of folds 


n A = 4.92 

[2 6 18] 

n 2 = 1.64 


n 3 = 0.55 


n\ = 2.97 

[51012] 

n' 2 = 1.48 


n' 2 = 1.24 


n\ = 2.99 

[15 30 36] 

ri 2 = 1.49 


n 2 = 1.25 


Table 11. Numbers of folds for baseline design. 


Determine the required antenna spacing, d,, for each channel. 


The required antenna spacing for all moduli sets are given below in 
Table 12. 


Moduli Set 

Antenna spacing 

[2 6 18] 

c / 1 =30.73 cm 
d 2 =10.24 cm 
c / 3 = 3.41 cm 

[51012] 

d\ =39.33 cm 
d 2 =19.67 cm 
c /3 =16.39 cm 

[15 30 36] 

d\ =39.33 cm 
d 2 =19.67 cm 
c /3 =16.39 cm 


Table 12. Antenna spacing for baseline design. 
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6 . 


Re-mapping the Field of View (FOV). 

The real moduli set’s FOV was not scaled and hence the scale 
factor % rea i =1. The FOV of the three-channel system is given in 
Table 13. 



Table 13. Scale factor and FOV for baseline design. 


7. Determine the thresholds. 

The threshold levels for all the moduli sets are shown graphically in 
Figure 25-27. 



m 3 =18 


Figure 25. Threshold levels for real moduli set [2 6 18]. 
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m' 3 = 36 

Figure 27. Threshold levels for second virtual moduli set [15 30 36]. 
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8 . 


Calculate the phase adjustment term, £ 7 , for each channel. 


The phase adjustment terms, 47 , for the moduli sets are listed in 
Table 14. 


Moduli set 

Phase Adjustment term, £ 


Cr 

= 2.36 

[2 6 18] 

<2 

= -1.48 


<3 

= 1.54 



= -0.94 

[51012] 

C 2 

= -2.15 


a 

= 1.79 


a 

= -1.15 

[15 30 36] 

& 

= -2.18 


a 

= 1.82 


Table 14. Phase adjustment term for baseline design. 


The ideal transfer function for the baseline design is given in Figure 28. 
The baseline design has a broadside resolution of 0.19 deg RMS. In practice, it 
is not possible to achieve this accuracy as there would be errors which have not 
been accounted for in the current baseline design. A one watt emitter at a 
detection range of 12 km for the current hardware is used as a scenario. This 
scenario gives a SNR of approximately 38 dB. Using this scenario and with the 
addition of Gaussian noise as described in Chapter VI, the transfer function is 
given in Figure 29. The RMS error is approximately 0.5 deg RMS. 
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Figure 28. Transfer function for the baseline design. 



Figure 29. Transfer function with noise. 


The transfer function shown in Figure 29 represents a snapshot of the 
transfer function in time. Over time, the transfer function would vary due to 
random noise. In order to figure out how the noise would affect the probability of 
estimating the correct AOA over time, a Monte Carlo simulation was carried out. 
For the Monte Carlo simulation, the following parameters were used: 
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Relative error, a = 3%. 


Required outcome probability , a = 90%. 

For the given scenario, there is a greater than 50% probability that the 
correct AOA will be estimated within the FOV as shown in Figure 30. In order to 
satisfy Equation (6.14), the number of times the simulation needs to be run is 
2777. In practice, 2777 readings would need to be taken from the emitter when 
the determining the AOA for the results of the Monte Carlo simulation to hold 
true. As this is not realistc, a comparison between the probabilities of taking a 
reduced sample size of 100 readings against 2777 readings is given in Figure 31. 
It shows that a reduced sample of 100 readings has good correlation to the full 
Monte Carlo simulation results. Therefore, even if a reduced sample size of 100 
readings is used, the results from the Monte Carlo simulation would still hold true. 



Figure 30. Probability of correct AOA estimation for 

baseline design. 
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Figure 31. Probability of correct AOA estimation for 
full sample size vs 100 sample size. 

B. EFFECTS OF SPACING ERROR 

As it is unlikely that the array would be precisely constructed, the effects of 
spacing error were investigated. By adding a manufacturing error of 2 mm and 
4 mm along the array axis into the element spacing, the transfer functions were 
plotted in Figure 32. The spacing error does not affect the transfer function 
significantly for spacing errors of 2 mm while spacing errors of 4 mm introduce 
more fading regions away from the broadside. These spacing errors correspond 
to 0.016 and 0.032 wavelength displacements in the horizontal direction. Note 
that position errors along the array axis do not affect AOA estimate near the 
broadside. Position errors out of the array plane will affect the AOA estimates 
near the broadside. 
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Figure 32. Transfer function with noise and spacing errors. 

The probability of correct AOA estimation for spacing errors of 2 mm and 
4 mm is given in Figure 33. Comparing Figure 30 and Figure 33, spacing errors 
do not affect the probability of correct AOA estimation. These results suggest 
that small spacing error does not have significant effect on the performance of 
the RSNS virtual spacing DF system. 



Figure 33. Probability of correct AOA estimation with spacing errors. 
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C. EFFECTS OF PHASE ERROR 

There may be residual phase errors introduced into the DF system that 
may not be possible to remove during the calibration. Such phase errors could 
result from a different antenna temperature and cable phase error. By adding a 
phase error of 5 degrees and 10 degrees to all the channels, the transfer function 
is plotted in Figure 34. Comparing Figure 29 with Figure 34 shows that the 
phase error affected the transfer function accuracy significantly. As the phase 
error increases, the transfer function has more fading regions. 



Figure 34. Transfer function with noise and phase errors. 

The probability of correct AOA estimation with phase errors is plotted in 
Figure 35. It showed with phase error the probability of correct AOA estimation 
reduces and increasing the phase error further reduces the probability. These 
results indicated that phase errors need to be minimized as much as possible 
during the hardware design as well as during DF system operation to improve the 
accuracy of the DF system. 
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Figure 35. Probability of correct AOA estimation with phase errors. 

D. EFFECTS OF FREQUENCY CHANGE 

In all the previous simulations, it was assumed that the frequency of the 
signal was fixed at 2.4 GHz. In order to investigate into how sensitive the 
baseline design was to changes in the frequency, the physical spacing of the 
antenna elements was fixed at ^ =30.93 cm, d 2 =10.44 cm and d 3 = 3.61 cm. 

The signal frequency was varied to investigate its effect on the performance of 
the baseline design. It was observed in Figure 36 that tolerance up to ±200 
MHz bandwidth is possible without affecting the accuracy of the transfer function. 
As the signal frequency deviates away from the design frequency, the 
corresponding transfer function also deviates accordingly. This can be seen in 
Figure 37 where the slope of the transfer function for a 2.0 GHz signal no longer 
resembles the slope of the transfer function for a 2.4 GHz signal. The change in 
slope is predictable and can be compensated in the processing. 
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Figure 36. Transfer function with noise using 2.4 GHz and 2.2 GHz signal. 



Figure 37. Transfer function with noise using 2.4 GHz and 2.0 GHz signal. 


As Figure 38 shows, there is no change in the probability of correct AOA 
estimation even with a 2.0 GHz signal. To design for a wideband frequency 
RSNS DF array would require several antenna baselines to cover the desired 
bandwidth with the necessary accuracy. 
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Figure 38. Probability of correct AOA estimation baseline design 

with 2.0 GHz signal. 

E. ALTERNATE MODULI SET 

The baseline design’s low value moduli set was selected after taking into 
consideration quantization errors in a low SNR scenario. The drawback of such 
a low moduli set was a small dynamic range, which resulted in lower resolution. 
To understand the effect of moduli change, an alternate higher moduli set was 
chosen for comparison while keeping the same h ratio. Table 15 shows the 
alternate moduli set and some of the key parameters obtained. 
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Moduli set 

Dynamic 

Range 

Antenna spacing 

FOV (degree) 

[4 12 36] 

M Real ~ 3 

c/-| =29.43 cm 
d 2 =9.81 cm 
c / 3 = 3.27 cm 

±90 

[10 20 24] 

M'virtuatt = 179 

d\ = 39.24 cm 
d 2 =19.62 cm 
c /3 =16.35 cm 

±28.4 

[20 40 48] 

MVirtuatZ = 359 

d\ = 39.24 cm 
d 2 =19.62 cm 
c /3 =16.35 cm 

±28.5 


Table 15. Alternate moduli set’s parameters. 


The alternate moduli set has an error of approximately 0.09 deg RMS 
against the previous baseline design error of 0.5 deg RMS. The increase in 
accuracy was due to the higher dynamic range of the new moduli set. Using a 
similar scenario, the transfer function and the probability was plotted in Figure 39 
and Figure 40. During multiple simulation runs, the occurrence of gross error 
was observed more often than the baseline design. This observation showed up 
in the resulting lower probability for correct AOA estimation in Figure 40. While 
the higher moduli set increases resolution, it trades off on the overall probability 
of correctly predicting the AOA. 
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Figure 39. Transfer function of alternate moduli set. 



Figure 40. Probability of correct AOA estimation for alternate moduli set. 

The effects of spacing errors and phase errors were also plotted for the 
alternate moduli set in Figure 41 and Figure 42. The tolerance to spacing error 
increased but the tolerance to phase error decreased. Using a 2.0 GHz signal, 
the alternate moduli set exhibits the similar effects of frequency changes as the 
baseline design. This can be seen in Figure 43. 
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Figure 41. Transfer function for alternate moduli with spacing errors. 



Figure 42. Transfer function for alternate moduli with phase errors. 
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Figure 43. Transfer function for alternate moduli using 2.4 GHz 

and 2.0 GHz signal. 

F. EFFECTS OF SNR CHANGE 

If the demodulator board is replaced with another of higher sensitivity, the 
baseline design would be able to handle a much lower SNR. This would enable 
the RSNS virtual DF system to increase its maximum range of detection. To 
investigate the effects of such low SNR, a one watt emitter is simulated at 30 km 
and 60 km. The probability of correct AOA estimation was plotted for both the 
baseline design and the alternate moduli set design. While both designs showed 
a reduction in probability, the alternate moduli set’s performance at longer range, 
which translates to low SNR, was inferior to the baseline design. This can be 
clearly seen by comparing Figure 44 and Figure 45 for a signal at 60 km. The 
probability of correct AOA estimation for the baseline design was more than 50% 
while it was less than 50% for the alternate moduli set. 

In order to select the correct moduli set, the most probable SNR condition 
needs to be assessed. If a low SNR condition is expected, the baseline design 
could be used. If a high SNR condition is expected, the alternate moduli set 
should be used instead. 
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Figure 44. Probability of correct AOA estimation for baseline design with 

emitter at various ranges. 



Figure 45. Probability of correct AOA estimation for alternate moduli with 

emitter at various ranges. 

G. USING TWO SET OF MODULI 

Since the RSNS-based interferometer was relatively insensitive to small 
spacing error, an investigation into using two complete sets of different moduli 
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based on the same ratio was carried out. Based on the ratio, h = ^, the 

baseline moduli set and alternate moduli set was selected. Using the baseline 
design’s spacing distance, the alternate moduli set was used together with the 
baseline moduli set. Although this effectively introduced spacing error into the 
alternate moduli set, the limited impact on the accuracy of the transfer function 
may be worth the tradeoff for higher probability of correct estimation in low SNR 
conditions. The transfer function for the baseline and alternate moduli set is 
shown in Figure 46. 



Figure 46. Transfer function for baseline and alternate modui set. 
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VIII. CONCLUSION 


This thesis focused on the design and analysis of a three-channel RSNS 
virtual spacing DF system. By analyzing a chosen design with addition of 
Gaussian noise as well as carrying out Monte Carlo simulations, the sensitivities 
to spacing error, phase error, frequency error, moduli set change and SNR 
change were determined. The baseline design showed good tolerance to most 
of these errors and should be considered for situations where low SNR is 
expected. For situations where high SNR is expected, the alternate moduli set 
should be selected. As the RSNS-based DF system was robust to the small 
spacing error, the best solution could be the dual use of both low and high moduli 
set. This would give a high probability of correct estimation at all time while 
retaining the accuracy of correcting mapping of the estimated AOA to the actual 
AOA. 

The application of RSNS virtual spacing to DF has been demonstrated 
both theoretically as well as experimentally in the past. It showed great promise 
in applying this concept to platforms that have limited size but yet require high 
accuracy for DF. From a military perspective, an unmanned aerial vehicle used 
to detect a low probability of intercept emitters is one of the likely scenarios 
where RSNS-based DF concept could be applied. 

Some of the possible future research areas to continue this RSNS-based 
DF concept are: 

1. To build the baseline design using available hardware and verify 
the performance of the design with its theoretical performance. With 
the experimental data, other errors could be identified and 
corrected. Initial testing could be done in the laboratory before 
expanding it to field testing. 

2. To determine the best set of two or more moduli set in order to 
optimize the performance of the RSNS-based DF system. 
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3. Further improvement of the simulation by including the effects of 
mutual coupling interference. 

4. The use of non-PRP moduli set resulted in a smaller dynamic range 
than what would have been possible with a PRP moduli set. 
Studies into the use of PRP moduli set for virtual spacing should 
continue in order to fully exploit the advantage RSNS-based DF 
system. 
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